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p Abstract: Robustness of the Luttinger theorem for fermionic hquids is examined 

in holography. The statement of the Luttinger theorem, the equahty between the 
fermion charge density and the volume enclosed by the Fermi surface, can be mapped 
psj to a Gauss's law in the gravity dual, a la Sachdev. We show that various deforma- 

QQ tions in the gravity dual, such as inclusion of magnetic fields, a parity-violating 

OO ^-term, dilatonic deformations, and higher-derivative corrections, do not violate the 

ly-^ holographic derivation of the Luttinger theorem, as long as the theory is in a con- 

(^ fining phase. Therefore a robustness of the theorem is found for strongly correlated 

fermions coupled with strongly coupled sectors which admit gravity duals. On the 
other hand, in the deconfined phase, we also show that the deficit appearing in the 
>■ Luttinger theorem is again universal. It measures a total deficit which measures the 

S^ charge of the deconfined ( "fractionalized" ) fermions, independent of the deformation 

^ parameters. 
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1. Luttinger theorem and holography 

The Luttinger theorem [1, 2] is one of the key fundamental relations in condensed 
matter physics and it states that the volume enclosed by the Fermi-surface is equal to 
the charge density. This theorem, which is originally derived by Luttinger and Ward 
for Landau's Fermi liquids, is non-trivial in the sense that this theorem is about the 
volume enclosed by the Fermi-surface. Remember that in the Landau's Fermi-liquid 
picture, we have a quasi-particle description for the spectrum near the Fermi-surface, 
but generically the quasi-particle description is not valid for the spectrum far away 
from the Fermi- surface, therefore the spectrum deep inside the Fermi-surface does 
not always allow the quasi-particle description generically. The non-trivial point of 
the Luttinger theorem is that it relates the spectrum not only near the Fermi-surface 
(where quasi-particle picture holds) but also deep inside the Fermi-surface (where 
quasi-particle picture does not hold), to the charge density. 

It is widely known that the theorem holds for Fermi liquids having a Fermi 
surface, and there is a general non-perturbative proof of the Luttinger theorem for 
Fermi liquids [3] (the original proof by Luttinger and Ward was with perturbation 
of Fermi liquids). The proof of [3] is based on a U{1) gauge symmetry, Fermi-liquids 
description near the Fermi-surface and a mild assumption for dynamical degrees of 



freedom, namely, all the momentum and charge carrying degrees of freedom are quasi- 
particles near the Fermi-surfaces. See also [4, 5] for further developments concerning 
the proof of the Luttinger theorem. 

On the other hand, in nature there are quite interestingly materials, such as 
high Tc superconductors or heavy fermions, where its normal phase shows non-Fermi 
liquid behavior, and in addition, the standard quasi-particle description breaks down. 
In such situations whether the theorem holds or not is to be better understood. 

Recent progress in applications of string theory, the holographic principle [6, 7, 
8], to condensed matter systems brought an insight about the Luttinger theorem 
in strongly correlated fermion systems. In [9, 10, 11], it was pointed out that a 
holographic system with a charged bulk fermion, which is called "electron star" in 
the literature, exhibits the Luttinger theorem of the boundary fermion theory. This 
is based on the observation that the bulk fermions obey a bulk Luttinger theorem 
at each radius for the electron star. Then, Hartnoll pointed out [12] that the flux 
emanating from the black hole horizon will equal the deviation from the Luttinger 
relation. Furthermore, Sachdev clarifled [13] that in a simple holographic set-up 
for fermions with fermion-number chemical potential, the Luttinger relation follows 
simply from the Gauss's law in the bulk and that it holds in conflned phase (thermal 
gas phase) but breaks down in deconflned phase (black hole phase). However it 
is also true that their argument uses a speciflc holographic setup like neglecting 
higher derivative corrections. Therefore it is natural to ask how universal the non- 
perturbative Luttinger theorem is for fermions. 

The holographic principle has been widely applied to various gravity setups, 
and robust correspondence has been thoroughly studied. Among many variations 
of the holographic models, some of the most popular and meaningful ones are: (i) 
higher-derivative corrections in the bulk gravity -|- Maxwell theory, (ii) inclusion of 9 
term and magnetic fleld, (iii) inclusion of a dilaton to have dilatonic gravity models.^ 
Each corresponds, in terms of condensed matter theory language, to: (i) Sub-leading 
terms concerning the strong coupling expansion, (ii) Parity-violating terms inducing 
quantum Hall effects under magnetic flelds, and (iii) Drastically different infra-red 
behavior, for example having a Lifshitz-like scaling near quantum critical points, and 
more realistic systems with vanishing entropy at zero temperature. 

We would like to study whether the holographic derivation of the Luttinger 
theorem a la Sachdev can survive against the deformations, to flnd a universality 
of the holographic Luttinger theorem. In this paper, we examine these popular 
deformations and show the Luttinger theorem to hold for all of these deformations, 
in the case of conflning phases.^ 

^For the inclusion of the dilaton, see also [14, 15, 16]. 

^The Luttinger theorem for various string-motivated field theories was studied in [17]. 



2. Holographic derivation of the Luttinger theorem 

We follow the beautiful argument of Sachdev's holographic derivation of the Luttinger 
relation [13], to show the robustness of the holographic Luttinger theorem. We 
generalize Sachidev's derivation [13], in particular concerning the following points: 
(i) In the gravity side we allow for a generic action for the bulk gauge fields including 
higher-derivative corrections, (ii) We allow an axion coupling (responsible for a 9- 
term) in the action and we include a magnetic field in the background, (iii) The 
gravity action is coupled to a dilaton with a generic form of its couplings to the 
gravity and the bulk gauge fields, resulting in generic background geometry which is 
different from the cut-off AdS4 used in [13]. Specifically, we consider the following 
action 

S = Jd'x^ [R - 2{V^f - g{^)Vo - {Vaf - /(</))£[Fm^] 

+ t^PiV^'DM + m)ij - h{<P) a e'^'^'^'^FKLFMN] ■ (2.1) 

Here, the spacetime dimension is 4 which is dual to a condensed matter system 
in 2-1-1 dimensions. The 2-1-1 dimensions are spanned by x,y and t, while the z 
direction is the emergent space coming out of the holographic principle. The sub- 
spacetime 2; = is the boundary of the bulk geometry. The geometry typically is an 
asymptotic AdS^ geometry, but in this paper we do not rely on any specific metric.^ 
D is the Dirac operator in which the f/(l) charge of the fermion is included as q in 
the covariant derivative. The Maxwell field in the bulk can have a generic nonlinear 
electrodynamics Lagrangian £.[F]. For example, the standard Maxwell Lagrangian is 
C[F] = (l/4)F/v,/7vF^^^, while the famous Dirac-Born-Infeld action which is natural 
in string theory is 



£{F) = ^ / det Lmn + \fmn] , (2.2) 

and includes higher derivative terms as a form of F^ and higher multiples. In hologra- 
phy, higher-derivative corrections may be related directly to a physical consequence; 
for example, the famous calculation on the shear viscosity of the quark-gluon plasma 
[18, 19] can be corrected by higher-derivative terms to have lower values, but a physi- 
cal constraint on the form of the higher-derivative terms coming from a bulk causality 
may forbid the value of the viscosity (divided by the entropy density) to go lower 
[20, 21]. In generic holographic setting, even if we take large A^ limit, these higher 
derivative corrections are non-negligible. 

The system couples to the bulk scalar field which is a dilaton in string theory. 
Explicit solutions (such as the ones with horizons studied in, for examples, [22] - [34] 



^The Luttinger theorem is in principle a low energy phenomena, so the UV behavior is expected 
not to be relevant for the discussion. 



) are not necessary in the following. In the second line of the action (2.1), we have 
the parity- violating ^-term (the axion coupling), which is relevant for a quantum 
Hall effect under the magnetic field, see for examples, [35, 36, 37, 29, 38]. 

The essence of the Sachdev's derivation is to consider a bulk fermion one-loop 
path-integration in the total free energy as a semi-classical approximation. This 
affects the scalar potential of the bulk gauge field, to minimize the free energy (which 
is nothing but the on-shell effective action in the gravity side). A radial integral of 
the bulk Gauss's law turns out to be nothing but the Luttinger relation. Once one 
puts (f) = a = and C = (1/4)F^, and takes a cut-off AdS space in the derivation 
below, everything reduces to the Sachdev's original derivation. 

The bulk free energy per a unit volume in our case is 

:F = jdz^ {-f{<P)C[F] - h{<P)aFF) - ^Tr Log [D-V + m]. (2.3) 

We have integrated out the bulk fermion '0, to have the last term.^ We ignore the 
back reaction of the fermions to the bulk geometry and the dilaton and the axion; we 
only consider a back reaction to the gauge field (a possible justification is presented 
in the next section). 

As for the gauge field configuration, we assume the homogeneity and the isotropy 
in the (x, y) directions. Then non-vanishing components of the static gauge field 
strengths are only Ftz and F^y Using the Jacobi identity in the {x,y,z) space, we 
obtain dzF^y = which means that B = F^y is constant. So, in the A^ = gauge, 
we are left with the gauge configuration 

Ao = Ao{z), A, = --y, Ay = -x, A, = 0. (2.4) 

The charge density of the boundary theory is given by 

{Q) = -^ = -Dziz = 0) where D, = f{<P)^^ + h{^)aF,,y. (2.5) 

Here we have defined the chemical potential as ;U = $(2; = 0) where we denote the 
temporal component of the gauge field as Aq = «$. Note that the $ dependence in 
the fermion loop in (2.3) does not contribute to the definition of the charge (2.5), 
because the bulk fermion wave functions vanish at the boundary z = Q where the 
chemical potential is defined, due to the normalizability. 

To explicitly perform the fermion one-loop integral to evaluate (2.3), one just 
needs a formal expression characterizing the discrete energy eigen modes of the bulk 
fermions, which looks 

T^n XLn{z) = Ei^n Xl,n{^)- (2-6) 



^If wc do not include this fermion path integral to account for the back-reaction to the gauge 
potential, we would not obtain the Luttinger relation (see for example a discussion in [39]). The 
effects of the bulk fermions should be communicated with the gauge potential, as the Luttinger 
theorem is a relation between the Fermi surface and the charge density. 



The operator P„ is nothing but a spatial part of the covariant Dirac operator in 
the curved geometry. xi,n and Ei^n are the eigenf unctions and eigenvalues of the 
operator D„, respectively.^ Here the integer / labels the Kaluza-Klein modes of the 
bulk fermion in the curved space, and n labels the Landau levels of the fermion wave 
function in the x-y space as the magnetic field is present. Below, we shall use only 
the two facts: First, the operator D„ includes a trivial term g$ coming from the 
minimal coupling in the Dirac operator, and second, the normalization is given as 



dzJ^^xU^)XiA^) = l, (2.7) 

V ~9tt 

where no summation for / and n is imposed.^ See Appendix A for the explicit 
evaluation of the states with the operator V. 

Using this energy eigenvalue, generically the trace log term in the free energy 
can be evaluated as 



TrLog [D-T + m] = ^Y,Y,Ei,n9{-EiA. {21 



The step function is necessary to count only the residues appearing in the shift of the 
poles in the path integral in the off-shell /c° space. See Appendix B for the detailed 
calculations. The factor qB /2tt is the unit volume of the discretized momentum 
space {kx, ky) due to the magnetic field. 

Now, to find a saddle point of the free energy with respect to the bulk field $, 
we consider the bulk on-shell equation 

To calculate this with the fermion loop term, we use (2.7) and (2.6) to rewrite the 
fermion free energy (2.8) as 

Ei^J{-EiA = jdzJ^^xU^)Xi,n{^)Ei,ne{-EiA 

= l^^'J?^, xU^)VxLn{^)0{-EiA . (2.10) 

With the fact that $ dependence of the operator V is just linear in g<l> as it is a 
Dirac operator (see Appendix A for the details), the minimization of the free energy 



^For deconfined geometry / can be a continuous parameter, but here we formally write the generic 
label as I. In addition to that, in the deconfined geometry, the energy eigenvalues are generically 
complex, thus the amplitude of the wave function damps in time exponentially, as in the case of 
quasi-normal modes. Here formally we regard our calculation performed in a time scale shorter 
than the decay time scale. 

^We consider the cases where bulk metric takes the form as ds^ — gtt{z)dt^ + X]i=i 9ii{z)d^i + 
g^^{z)dz'^. 



(2.9) is^ 

-5.(/(0)^ + M0)«i^..) -?|fy^E^(-^'.X«(^)^^ (2.11) 

Then we make an integration over the z space, which leads us to the Luttinger 
relation, as follows. For that, we need one more information for the geometry at 
the IR. Let us first consider a generic confining geometry, for which the geometry 
consistently ends at 2; = zy^. The Gauss's law at the IR end z = zi^i shows that the 
electric flux Dz vanishes there. Then, from the z integration of (2.11), we obtain 

(Q)/l=^^J2^i-Ei,n). (2.12) 

l,n 

This is the Luttinger relation, since the right hand side is the volume enclosed by 
the Fermi surface. Note that since we turn on the magnetic field, the Landau levels 
appear and the x-y momentum space is discretized, and resultantly the unit volume 
of the x-y momentum space qB/2'K appears. 

For deconfined geometries, the IR boundary condition differs, and in particular 
the electric flux does not vanish there. There appears a deficit in the Luttinger 
relation. We discuss the situation in the next section. 

One should have noticed that the derivation here is almost identical with what 
Sachdev gave in [13]. However, we find it intriguing that the derivation by Sachdev 
is so robust that the theorem is valid against various deformations of the system. 
In particular, the inclusion of the higher derivative corrections corresponds to the 
direction toward a weak coupling where it is plausible that the Luttinger theorem is 
valid. Furthermore, the inclusion of the background magnetic field is interesting, as it 
not only introduces a nice regularization of the momentum space but also is involved 
with quantum Hall effects. The dilatonic corrections are related with existence of 
different scaling at the IR, and even in those systems our generalized derivation shows 
that the Luttinger theorem holds. 

Although in this paper we worked in four spacetime dimensions in the grav- 
ity theory (which corresponds to three spacetime dimensions for the liquid system), 
we can generalize the derivation to higher dimensions, in a straightforward man- 
ner. A possible obstacle would be the F A F term in the bulk, which should be 
generalized to a Chern-Simons term in higher dimensions, then one cannot impose 
the isotropy which we have employed in the derivation above. Another concern 
may be on the fermion integral, since generically in higher dimensions the theory 
becomes non-renormalizable. However, the leading quantum loop which we consid- 
ered in this paper is just an effect of the chemical potential of the one-loop diagram 

^When taking a variation of E9{—E) with respect to $, one may be worried about the $- 
dependence in E inside the step function. However, as the variation of 0{E) is a delta function 
while there is an overall E in front of it, the worring contribution disappears. 



of fermions, which can be defined without any problem. Therefore we claim that 
higher-dimensional generalization of the derivation can be done accordingly. 

3. Discussions 

The essential statement which we would like to make in this paper is just the robust- 
ness of the holographic Luttinger theorem, which was already shown in the previous 
section. We end this short paper with two important observations: First, the validity 
of our calculation in view of quantum corrections in string theory and AdS/CFT, 
and second, the emergent dependence on deformation parameters only in the decon- 
finement phase. 

3.1 Quantum corrections in string theory 

In the derivation of the holographic Luttinger theorem, Sachdev and we computed 
the fermion one-loop diagram and considered its back-reaction to the gauge potential. 
However, in general in string theory, other loops involving gravity and other fields 
may contribute, so the effect on the gauge potential may not be only from the fermion 
one-loop. Here we shall point out that a proper large N scaling in AdS/CFT can 
avoid this mixing problem*. 

We assume that the bulk fermion is from a space-filling D-brane. This fermion 
is often called "mesino" since it may be a fermion counterpart of mesons, when the 
D-brane is identified with the fiavor D-branes in the AdS/CFT correspondence [40]. 
The gauge potential A^ also comes from a space-filling D-brane. This means that the 
effective gauge coupling q is of order 0{1/\/N), because in the AdS/CFT correspon- 
dence, the string coupling constant Qs scales as ~ l/N, and the coupling on the D- 
brane is the open string coupling y/g^ while the coupling in the bulk geometry (gravity 
and the dilaton (p and the axion a) is Qs- Denoting the graviton/dilaton/axion fluc- 
tuation as 6g, and the gauge fluctuation as SA^, and the fermion fluctuation as ip, 
then the generic dependence in N in general AdS/CFT is written as 

S = {d6gr + j^{69f + ^^{69)' 
+ {dSAM? + ^T^Dm^ + -^(Mm)# + ^<^^(Mm)' + ^<5#^ + ■ ■ ■ (3.1) 



Looking back our derivation of the holographic Luttinger theorem, we have performed 
the if) integral in this action. From the generic action written above, we observe that 
this one-loop integral involves the term -t={5A^)iI)iIj, so it shifts the action by a 
term of order 0{\/yN). This is a leading order effect, compared to the quantum 



^See also [13] for the argument to include quantum fluctuations beyond the fermion one-loop. 



corrections involving the graviton, the dilaton and the axion,^ since in the action 
above those corrections come with a couphng of 0{1/N). 

Therefore, we conclude that, assuming that the origin of the fermions and the 
gauge fields is D-branes in the bulk, the fermion one-loop integral is the leading order 
in 1/A^ expansion in AdS/CFT, so we can consistently ignore the other quantum 
corrections. 

3.2 Deconfined phase 

As first pointed out by HartnoU [12] and also studied in detail by Iqbal and Liu 
recently [15], the presence of the black hole horizon significantly alters the result; in 
the final line of the derivation of the Luttinger theorem (2.12), we have used the fact 
that the electric displacement D^ at the IR endpoint vanishes due to the confining 
geometry. However in the presence of the black hole the IR boundary condition is 
different and there exists in general an electric flux emanating from the black hole 
horizon, for charged black holes. 



(Q) + ^.|_,^^=gE^(-^'--)' (3.2) 



Ln 



with 

So, there exists a Luttinger deficit for the deconfined phase [12, 15]. Here D^ is the 
electric displacement in the bulk at the IR of the geometry, which is the electric flux 
penetrating the horizon of the black hole. This comes from the first term of (2.11). 

In the presence of the deformations which we consider, for the deconfined phase, 
there appears a dependence on the magnetic field and the axion field at the horizon. 
The axion field corresponds to the parity-violating 9 term, so in general, the Luttinger 
theorem in the deconfinement phase is violated with a deficit dependent on the 9 term 
and the magnetic field in addition to the deficit [15] of the electric charge carried 
by the "fractionalized" fermions [41, 42, 17] (which are deconfined quarks in the 
standard holographic QCD terminology). 

The electric displacement includes a contribution from the parity- violating term, 
the second term of the (3.3), 



D. 



= 9,sB. (3.4) 

parity— odd 



^Here, as a classical background geometry, we have assumed that the back-reaction of the flavor 
brane itself (with t/i = 0) is already included in the geometry. In this paper we need not to specify 
the geometry for the derivation of the Luttinger theorem, so whatever the back-reacted classical 
geometry is, there is no problem in the derivation. 



This term directly responds to the magnetic field present in the system, and the 
effective value of the coefficient, 6es, in the Luttinger theorem is given by 

des = [^(0)a]2=2iR- (3.5) 

It is intriguing that the parameter 6es is not given by the UV geometry but the IR 
geometry. The IR geometry is not directly related to the parameters of the fermion 
liquid system defined at UV, and it is rather determined by the gravitational dynam- 
ics corresponding to strongly coupled sectors, which these fractionalized fermions 
couple. This dependence on IR geometry refiects the fact that the Luttinger relation 
is a phenomena at low energy. Note that even with no electric field Fqz, once the 
parity-violating term and the magnetic field F^y is turned on, the Luttinger deficit 
appears. 

The deficit appearing in the Luttinger relation (3.2) appears to depend explicitly 
on the deformation parameters we introduced in the gravity dual. The parameters 
are in the nonlinear electrodynamics C[F] and the dilaton-axion couplings to the 
gauge fields in the gravity side. However, in (3.2), the deficit depend only on the 
electric displacement D^ at the black hole horizon. The effect of the parity-odd 
term {9csB) is also included in the electric displacement. So, we conclude that the 
effect of the deconfinement phase to the deviation from Luttinger theorem can be 
summarized into electric displacement Dz, which is determined by the total charge 
of the "fractionalized" fermions, even in the presence of the deformations considered 
in this paper. Note that we have not assumed the existence of the quasi-particle 
picture in our derivation. 
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A. Bulk fermion eigen-states 

In this appendix we shall explicitly calculate (2.6) with the background gauge field. 
We work in the Lorentzian signature, while the Euclidean signature (which we em- 
ployed in the derivation) can be easily obtained by an analytic continuation. First, 
we derive the Dirac operator. The fermion action in the bulk is 

^fermion = fd^^'x ^^ I [^r^'l^Af^ - 771^,^] . (A.l) 



Here the Dirac operator is Dm = Om + \wahM^°''^ — iq^M- The definition of the 
Gamma matrices in the local Lorentz frame are 

with 7° = icTs, 7^ = CTi, and 7^ = — cr2, where cri,cr2,0"3 are the Pauli matrices. We 
follow the notation of [43] and [32] except for the assignment of 7^^ (this difference 
is necessary to see the diagonalization as for fermion components, see below). The 
notation for the indices are: M = 0, 1, 2, 2, and /i = 0, 1, 2. 
Writing the 4-component fermion as 



^ = ( ^" ) , ^± = {-9gn~">±. (A.3) 




where (l)± is a two-spinor, the Dirac equation is 



{d, T m^^,)(f)± = TiK^-i^(p^, (A.4) 

9zz 

with A'o = —u/^^^{do — iqAo) and K"j = —i{di — iqAi) with z = 1, 2. 

In the Dirac equation, the four spinor components are coupled, while we would 
like to group them into 2-spinors to derive (2.6). In our convention, among the 
gamma matrices, 7° is a diagonal matrix while 7* is not, so if we can bring the -R'i7* 
to a diagonal form, the decomposition to the 2-spinors is complete. This is nothing 
but solving the following eigen equation in the x-y space, 

z(iri7^ + i^27')0±=(^|^°^<^±, (A.5) 

where a and /3 are complex constants, and </'±(x, y) are 2-spinor wave functions rep- 
resenting Landau levels of fermions in magnetic fields. The gauge field configuration 
(2.4) satisfies [Ki,K2\ = iqB, so using a creation and an annihilation operator 

{Ki + iK2) = a, -^=={Ki - iK2) = a\ (A.6) 



V'2^ ' 72^ 

we have the harmonic osccillator [a, a)] = 1. The eigen equation (A.5) becomes 

a\ ~ (a 



iv^2^( 



Vato;^^=Vo/3,'''^- ^"^-^^ 



A solution can be easily found as 



[zV2^57^/a)(at)"-i|0) 
:at)"|0) 



-^^(,.v^.^;™«; '^M, n = 0,l,2,.... (A. 



10 



We need /3 = —2qBn/a to have this solution. With this wave function for {x,y), we 
choose 



^±/ 



.^/.,„We— I ''^^/"t,r >"'">), (A.9) 



(at)''|0) 



then the Dirac equation (A. 4) is written as 



9. 



{dz =F m^Qzz) <P±,n 



zz 



From the phase matching in (A. 10), we have Wj^ = W- = En- The upper spinor 
component of the equation (A. 10) yields 



[Oz - m^/g^z) = - [En + qAo + J a_ , (A.llJ 



Qzz ot+ V V fi'jJ / " 



{Oz + m^/gTz) =\En + qAQ + J «+ . (A.12) 



Qzz a- \ V 9ii J tt+ 

The lower spinor component of the equation (A. 10) yields 



^" [dz - m^^z) f+A^) = - (-{En + gAo) + J-^(-2g5n/a_) ) /_,„(^) 

(A.13) 



^" {dz + m^/^z) f-,n{^) = ( -{En + g^o) + J-^{-2qBn/a+) ) /+,„(^) . 



fi'zz \ V 9i' 

(A. 14) 

For (A.13) and (A. 14) to be consistent with (A. 11) and (A.12), we need 

a+ = — a_, a_ = 2qBn/a-., a+ = 2qBn/a+. (A. 15) 

This can be solved as 

a± = ±^2qBn or a± = tV^qB^ ■ (A.16) 

Then (A.13) and (A. 14) become the same as (A. 11) and (A.12), and we obtain two 
independent solutions as (A.16). 

The resultant equations (A.13) and (A. 14) can be re-written as 

Vt'^Xn{z)=EnXn{z), (A.17) 



9ii 



11 



with X = {f+if-Y'- This is for the choice a± = ±y/2qBn. For the other choice 
a± = TV^qBn, 



jy{i=2) 



-9tt 



{-1028^ - oim^/g^^) - qAo + ^3 



-gtt 



2qBn. (A.19) 



9zz V 9ii 

The operator P„ of (2.6) appears here as Vn . See also [13, 14]. Note that due 
to the interaction term between spins and magnetic fields, there are two V^s labeled 
by i. Generically this eigen-equation allows only a discrete spectrum^*^ for En, so the 
energy En is also labeled by / (the Kaluza-Klein modes in z space) in addition to n 
(the Landau levels) and i (spins) as 



T>n Xn,l,i{z) = En,l,i Xn,l,i{z) ■ 



(A.20) 



For simplicity of the notation, we will omit the index i in this paper. 

The normalizability condition should single out a certain linear combination of 
the two solutions. The reason why we got two solutions generically can be easily un- 
derstood by the following argument. The Dirac equation (A. 4) is a coupled equation 
of 0+ and 0_, but one can eliminate one of them. Bringing (A. 4) into the following 
form formally. 



V 9zz 
V 9zz 

we can combine these to eliminate 0_, to have 



(A.21) 
(A.22) 



[iK^Y 





"iVfi'zzJ 



..(A.23) 



9ii /o I / \/ ■ T^ u\ — l 9i' 

— {o^ + m^/gZ)[-iK^Y) \ — 
9zz V 9zz 

This is a 2nd order differential equation for a two-spinor 0_, so we generically have 
four independent solutions. Without magnetic field, given E, four states are de- 
generate and they correspond to spin {up and down}, and {normalizable and non- 
normalizable} modes. In the presence of magnetic fields, this degeneracy splits up 
by spins, so given spin and energy eigenvalue E, there is a set of a normalizable and 
a non-normalizable mode. 



B. Evaluation of the fermion free energy 

In this appendix, we demonstrate the calculation of the fermion one-loop free energy 
(2.8). 

lOpor given index n and z, the corresponding parts of Dirac equation be coupled 1st order dif- 
ferential equations for two spinor components, so it allows two independent solutions. Taking the 
appropriate ration between these two, we can eliminate non-normalizable mode in the bulk UV and 
regular mode in bulk IR in Euclidian signature. Then it allows only discrete modes labeled by I. 



12 



B.l Reduction to 1 + dimension 

In the previous appendix, we demonstrate that given / (the Kaluza-Klein modes in 
z space), n (the Landau levels) and % (spins), the Dirac equation reduces to a single 
equation (A. 20) which determines -E/,^- Therefore the action for the fermion should 
read as 



S = J2fdt t^i^rrit) {tdt - Ei^n) ^«,n(t). (B.r 



l,n 

Here ^!in{t) is a single component fermion, since the dependence on the Landau level 
n is already included with the specific Landau level wave function in the x-y space. 
The fermion field '^i,n{t) is properly normalized to have the action above. Note also 
that there are two Ei^^s, depending on the choice of the Dirac operator Vn in (A. 18) 
and (A. 19) labeled by i. In the following, we omit the index i for simplicity of the 
notation. 

Let us proceed to calculate the free energy. We bring the Tr Log into the following 
expression, 

•^fennion - E |f E ^^ Log ^ - E,^) = ^ E ^ ^^ ^'W^) 

n I l,n "^-^^'i ^ ^ > 

= «^Vrd,/^^. (B.2) 



Ln -"^l 



Here, we have used the fact that the momentum integration in the kx-ky space is 
now replaced by the Landau level summation, as 

/dk^dky r 27rkdk f 7rd{k^) ^-^ 2'KqB ^^ qB 
-J2^ = jj2^ = j^2^ = 2^j2^ = 2^2^- ^ -3) 

We have used the momentum relation k"^ = 2qBn which we obtained in the last 
appendix. 

Now we perform the off-shell w integration. The standard path in the complex 
w-plane for the integration of w rounds the upper half plane plus the real axis. In the 
ie prescription, the pole contributing in the path integral is the one on the negative 
real axis of w. This appears only for the negative s so the integral is non-zero only 
when Ein < 0, and we obtain 

-^fermion = |f E / ds^27li9{-Ei,n) = -l^ E ^«,n^(-^«,n). (B.4) 

l,n -^^l," l,n 

This expression is (2.8) in the Euclidean notation. 
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B.2 Another viewpoint: Reduction to 1 + 2 dimension 

We shall present another viewpoint here, to evaluate the fermion free energy, which 
is a natural dimensional reduction along z. 

First, we shall decompose the bulk fermion into Kaluza-Klein modes along z 
labeled by /. The effective action for the /-th mode has a kinetic operator -D^7^ — vni 
where m; is the mass for the /-th mode of the decomposed fermion. Then, we compute 
Tr Log(D^7'^ — vdi) for each state / and make a summation over /. 

Let us work out the Kaluza-Klein decomposition explicitly. The bulk fermion 
action is 

^fermion = U'-^^X ^^ I [ ^F^^Dm^ - m^'V' ] • (B.5) 

Using the notation in the previous appendix, this action S'fermion can be explicitly 
written as 



£^^xJ^i (0V7°,0^-7° 




(B.6) 



Here 



K^ ^ J-^K^ = {-t{do - tqAo) , -tJ-^{d, - zqA,))) , (B.7) 

V 9ii V 9ii 

and Dz = dz — iqA^. In view of this, we consider the following matrix equation so 
that the above matrix is diagonalized. 



niAY) (B.8 




9zz 



where nii is some eigenvalue. We write^^ 

where \E' is a two-spinor which is a function of (t, x, y), and g±{z) are scalar functions. 
These g±{z) are required to satisfy 

yf^tm j \9-{z)) \9l'{z)J 

"'^^The reason why we took a common factor g±{z) for the 2-spinor (j)± is that this z-dependent 
factor go through the 7 matrices in K^^^^ such that the z integration can be done independently 
as (B.ll), (B.12). 
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(B.13) 



With explicit eigen wave functions, we can reduce the action to a 3-diniensional 
action. Substituting (B.9) to the action, we obtain 

l,m J \ 9tt 

(B.ll) 
The z integration gives a normahzation for the fermion which we require as 

dz^^[gf*{z)g^-\z)+g'l^*{z)g'^-\z))=6,m. (B.12) 

Furthermore, we calculate the chemical potential for each mode as 

f^im ^ jdz ^Mz) {gT{z)gt\z) + g^^*{z)g^-\z)) . 

Here we assume that this chemical potential is diagonal 

and also consider only the case below where gn{z) = —gu{z) such that dimensional 
reduction along z ensures the Lorentz invariance in the t-x-y spacetime^^. Then we 
finally obtain the decomposition 

^fermion = Yl J^' ^ ' ^^'\D .l" " mi)¥^ , (B.15) 

I ^ 

where Di = di — iqAi and Dq = do — iqfii- 

The on-shell condition for the l-th. fermion is Ei + g/i/ = ym|+X^ where k is 
the magnitude of the momentum in the {x,y) space. ^^ The derivative in the Dirac 
operator has eigenvalues d^ = {—i{w + qfii),iki, ik2). In the presence of the magnetic 
field, the momentum k is replaced by the Landau levels. 

Let us proceed to calculate the free energy. We bring the Tr Log into the following 
expression, 

•^fermion = ^^ Tr Log(D^7^ - m,) = ^^ Tr / dt Tr— — 

= E -fermion- (B-IG) 

I 

^^In general such diagonalization makes the diagonal matrix mi off-diagonal. So choosing the 
basis such that both /^;„, and mi diagonal is impossible for a generic Aq(z). Similarly choosing 
9tt{z.) = —gii{z) is impossible for a generic metric g^i, such as the Lifshitz form. The argument in 
the previous subsection using the 1 + 0-dimensional picture does not refer to these assumptions. 

^^This El coincides with Ei^n given in (A. 20), because both £"s, which is the eigenvalue of idt, 
are obtained from the equations of motion of same action; One is from 4-dimensional viewpoint 
(B.5), and the other is from dimensionally reduced 3-dimensional viewpoint (B.15). 
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As you can easily check, this integral is divergent. Therefore, we need to subtract the 
vacuum contribution (the free energy with ^i = 0) which corresponds to the Dirac 
fermi sea. Then 

-^ fermion = 'M'ermion (A'-U ^ -M'crmion lA''' — U) • [^■^') 

Now we perform the off-shell w integration. The standard path in the complex w- 
plane for the integration of w, which rounds the upper half plane, concerns two poles 
at u; = —qui ± y/k^ + f^. In the ie prescription, the only relevant pole in the path 
integral is the one on the negative real axis of w. However, note that now we have 
a contribution from the chemical potential, thus there is the case when the both of 
the two poles are on the negative real axis. 



/■oo 
J mi 



nil 

oo 

dt 

nil 



dwd'^k 2t 



d^k \ 2-ni 2t ^ , r— ^, 

(B.18) 



2y/W+¥ ^^^ 



Note that the 9 function in the first term is always equal to the unity. So, the first 
term is independent of /i^. The second term vanishes when /x; = 0. Therefore, we 
find that v^fermion — Ylii •^fcrmion coiucidcs with the sum of the second term in (B.18), 

t-^ /" d'^k /"°° 2f 

^feZion = J2^'J J^sj^dt 2vr^^^^==^(g/iz - VWTT'). (B.19) 

The 6 function represents the Fermi surface, since one satisfies the momentum con- 
straint relevant to the chemical potential, q^i > y/k"^ + t^ (roughly speaking, t is the 
mass of the l-th mode of the fermion). In fact, this is a widely known technique [44]. 
Denoting the value of t satisfying q^i — \/k^ + 1^ = as t*;, (B.19) is integrated to 



^-^ f a k C 2t 

= -' E / (0«<(*)9(-Bi(*)). (B.20) 



In the last equality, we have used the fact that t^i > mi is equivalent to y^k"^ + mf 
qfii < 0. 
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In the presence of the magnetic field, the momentum k'^ is replaced by the Landau 
levels 2qBn with a non-negative integer n. The momentum integral is accordingly 
normalized as dkidk2 = 2'n'kdk = TTd{k'^) = 27TqBdn, so 



I n 

This expression is (2.8) in the Euclidean notation. 
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